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Abstract. The evolution of atomic solitary waves in Bose-Einstein condensate (BEC) 
under adiabatic changes of the atomic scattering length is investigated. The variations 
of amplitude, width, and velocity of soliton are found for both spatial and time 
adiabatic variations. The possibility to use these variations to compress solitons up to 
very high local matter densities is shown both in absence and in presence of a parabolic 
C^) • confining potential. 
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1. Introduction 

The discovery of Bose-Einstein condenstates (BEC) in atomic vapors of alkali metals 
has opened new possibilities to observe matter waves of soliton type. These excitations 
were first discovered in BEC with repulsive interatomic interactions, in the form of dark 
solitons pQ, and more recently, in BEC with attractive interatomic interactions, in the 
form of bright solitons moving on zero backgrounds [2101 ■ Since soliton matter waves 
are of primary importance for developing concrete applications, it is of interest to devise 
methods which allow to control them. One possibility is to vary the atomic scattering 
length by means of external magnetic fields, i.e. by using Feshbach resonances. In 
Ref . [3] , it was shown that an abrupt change in time of the atomic scattering length, can 
lead to the splitting of the soliton with the generation of new solitons. The fragmentation 
of the wave obviously decreases the number of atoms contained in the original pulse, 
this being undesirable for applications such as atom lasers. 

In the present paper we suggest to use adiabatic variations of the atomic scattering 
length, both in time and in space, as an effective way for controlling soliton's parameters 
and to induce changes in their shape which could be useful for applications. In contrast 
to abrupt variations, adiabatic changes make possible to preserve the integrity of the 
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soliton (no splitting occurs), this leading for bright solitons to the compression of the 
pulse with the increase of the matter density, and for dark solitons to the compression 
of the hole with a decrease of the background level (matter from background move 
inside the hole). These phenomena are shown to exist both in presence and in absence 
of a parabolic confining potential. In the last case we derive analytical equations for 
soliton parameters in terms of adiabatic invariants and soliton perturbation theory. 
The results of this analysis is found in good agreement with direct integration of the 
ID Gross-Pitaevskii equation (GPE). The effect of a confining parabolic potential on 
the adiabatic dynamics of the soliton is also investigated by numerical simulations. We 
find that, for bright solitons, except for the oscillatory motion around the bottom of 
the trap, the phenomena of pulse compression is practically the same as in absence of 
the trap (this is particularly true for solitons initially at rest in the bottom of the trap). 
For dark soliton the phenomenon is only qualitatively preserved, due to the boundaries 
effects introduced by the trap. The possibility to compress BEC solitons could be an 
experimental tool to investigate the range of validity of the ID GPE. Since the quasi 
one dimensional regime is valid for low densities, it would be indeed interesting to see 
how far one can compress a soliton in a real experiment by means of adiabatic changes 
of the scattering length. Effects of adiabatic perturbations on soliton dynamics were 
also investigated in the context of Josephson junctions 5J and in nonlinear optics jE]. In 
contrast with Josephson junctions and optical fibers, which require structural changes 
or preparation of new samples, the study of adiabatic nonlinear perturbations on BEC 
solitons appears more natural and easy to perform, since the strength of the nonlinear 
interaction can be changed by using only external fields. We also remark that soliton 
dynamics in a quasi one-dimensional BEC under time- dependent linear potential was 
recently studied in [7j. 

The paper is organized as follows. In Section 1 we perform an analytical study of 
the effects on bright and dark BEC solitons induced by spatial and temporal adiabatic 
variations of the scattering length. For temporal variations we use both a variational 
approach and perturbation theory to obtain equations for soliton's parameters both in 
presence and absence of a confining parabolic potential. In the former case the problem 
can be linked to an adiabatic invariant of the Kepler problem with results depending only 
on adiabaticity (smoothness) but not on the smallness of the perturbation. For spatial 
variations we use standard soliton's perturbation theory. In Section 2 we compare the 
results of this analysis with direct numerical integrations of the full ID GPE both in 
absence and in presence of an harmonic trap. In the last section we summarize the main 
results of the paper. 
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2. Analysis 

The dynamics of a dilute trapped BEC is described by the Gross-Pitaevskii equation 
(GPE) 

ifuh{r,t) = -^A<j>(r,t) + VUr,t)(f>(r,t)+m 2 <f>{r,t), (1) 
2m 

with T = 47r ^ as , m being the atomic mass, a s the scattering length, and V tr = ^m(uj 2 x 2 + 
uj 2 y 2 + 003Z 2 ) represents the harmonic trap. The quasi ID geometry corresponds to the 
case {J 1 << cj| 3 . To model adiabatic variation of the atomic scattering in ID (cigar 
shaped) BEC we consider the following normalized Gross-Pitaevskii equation [315] 

iipt + *Pxx + <77(ar, t)V#| 2 - 0J 2 x 2 ip = 0, (2) 

where ip is the ground state wavefunction of the condensate, 7(2, £) is a slowly varying 
function of space and time, a ± 1 corresponds to the case of negative and positive 
scattering length a s {ui denotes the longitudinal frequency of the trap). In Eq. (J2J) the 
space has been normalized with respect to the healing length £ = 1/ ^8irp\a s \ (p is the 
atomic density), while the time with respect to to — rnt^/h. We remark that the ID 
approximation is valid for the number of atoms N < l T /\a s \, where l r = ^Jh/ (77^2,3) is 
the transverse oscillator length. It is also worth to mention that in the strong interaction 
limit a s N\ip\ 2 » 1, (a s > 0), the effective nonlinearity is of the type [E 3 so that 
deviations from Eq. (0) occur (this will not be considered in this paper). Although the 
analysis can be performed for a generic smooth functions j(x, t), we shall restrict to the 
limiting cases: 7 = 7(t), and 7 = 7(2), the former being experimentally more easy to 
realize. 



2.1. Case 7 = j(t) 

To study adiabatic temporal variations of the scattering length, case 7 = j(t), we 
remark that, due to the number of atoms conservation, it is natural to look for solution 
of Eq. (j2J) with a = +1 (negative scattering length) in the form of a bright soliton with 
time dependent parameters 

Mx,t) = A(t)sech( X ~ C , (t) )e^^-^ f )) 2+ ^^W) + ^, (3) 

a(t) 

where A(t), a(t), 4>(t), ((t), C(t), b(t), denote, respectively, amplitude, width, phase, 
center, velocity and " chirp" oscillation of the perturbed soliton. Taking these parameters 
(with their derivatives), as collective coordinates for the soliton, one can derive their 
time evolution from the Euler-Lagrange equations with respect to the space averaged 
Lagrangian L = J L(x,t)dx of Eq. (0) [10J. This gives for the center of mass: 
Qt + 4:LU 2 ( = 0, and for the width [TTj : 

2 16 4 7 (t)jV 
a tt + Au a= 5-5—. (4) 

Notice that Eq. (j3J) has the form of a perturbed Kepler problem ^2] • This analogy can 
be used to reformulate the problem in terms of the evolution of a unit mass moving in an 
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adiabatically varying Kepler potential. Since action-angle variables for the unperturbed 
case are known, one can describe the adiabatic change of the soliton width by means of 
the adiabatic invariant of the Kepler problem. As the result we find for uj 2 < 1 that 

4 64uj 2 

r l-7-77^d- (5) 



7 (t)iV l ( 7 (t)iV) 4 

We see that for uj = one recovers the law a = 4/( 7 (i)iV) known from nonlinear fiber 
optics. In this case, one can also show that the frequency of the width oscillation is given 
by uj a ~ 7 2 (t)iV 2 /47r. The analysis can be easily extended to two dimensions (2D) for the 
radially symmetric GPE. In this case the variational equation for the width becomes 
a tt + 2uj 2 a = 2(2 ~^ 7l(f)) , where Ni = N/{2-k)^ x = i/{y/2nl x ). Applying the above 
arguments, one finds that the width of a 2D soliton changes as a m (2 — N\^\{t)p* /u) x ' 2 . 
Eq. (jSJ) predicts for small uj a dependence I/7 of the soliton width . This result can 
also be obtained by perturbation theory. To this end we assume in the following that 
the range of variation of the potential is very large in comparison with the size of the 
condensate, so that we can neglect the uj term in Eq. ((2j) (this is particularly appropriate 
for static bright solitons located at the bottom of the trap, as we will show later). In 
this case the transformation ip = u/y/j (for 7 > 0), allows to rewrite Eq. (J2J) in the 
form of a perturbed nonlinear Schrodinger equation 

T 2 L 2 

iu t + u xx + \u\ 2 u = i(ln y/j) t u + (ln( 7 )) x .w x . + 0(^, -§) = R(u), (6) 

where T and L are, respectively, the temporal and spatial characteristic scales of the 
inhomogeneity, while T s , L s are characteristic time and space soliton scales, assumed 
to be much smaller than the corresponding inhomogeneity scales. In the case 7^ = 0, 
the right hand side of Eq. (JBJ) represents a small linear time dependent amplification 
(damping). If we consider the initial condition as a single soliton 

u(x,t) = V2A{t)sech{A{t){x - ())e ie{t) , 9 = k(t) + C(t)(x - C(t)), (7) 

then, using the equation for the energy N = J \u\ 2 dx, we find A = A j(t) this giving, 

for the amplitude and width of the soliton, A^ = A 0a /7, a = l/(A j), in agreement 

with the variational approach. From this we conclude that an adiabatic increase of 

the scattering length can be used to narrowing the width of a bright soliton matter 

wave. Soliton compression phenomena induced by linear damping amplification are 

also known from nonlinear optics [HIE]- A similar analysis can be performed for the 

adiabatic evolution of a dark soliton (a = — 1 in Eq.(J2J). In this case, the adiabatic 

evolution of the dark soliton ipd and of the background ub can be found as 

it I — /■< 
il) d = u tanh( — -^(x))e~ ld , u B = M 0v / 7exp(-z6'), 9 = -iu 2 Q / >y(t f )dt',( 
V2 Jo 

where Uq = y/JI, with \i the chemical potential of the condensate (notice that this last 

expression for ub is exact). From this expression we see that the width of the dark soliton 

is changing according to = a^o/ a/ 7 (£), i.e. a factor 1/ a/t(^) l ess m comparison with 

bright solitons. 



Adiabatic compression of BEC solitons 



5 




100; 




100 



Figure 1. a) Soliton's width vs. time for different values of A s . Panel a) refers to 
the bright soliton case (a = 1) The values of A s are stepped by 2, from A s = 2 (top 
curve) to A s = 10 (bottom curve). The other parameters are s = 0.5, Tf — 100, a = 1, 
v s = 0. An enlargement of the broken line box, showing rapidly decreasing oscillations 
for increasing A s values, is reported in the inset. Panel b) refers to the dark soliton 
case (er = —1 for A s values 2,4,8, (from top to bottom). The other parameters are 
s = 0.1, T f = 100, v s = 0. 



2.2. Case 7 = 7(2) 

The perturbative approach can be used to study also the case of spatial variations of the 
scattering length i.e. 7 = j(x). We have from Eq. (JHJ): R(u) = F(x)u x , F = (\n'j(x)) x , 
with R(u) considered to be a small perturbation. Using the perturbation theory for 
solitons [Hj, we find that the equations for the soliton's parameters (JJJ) are: 

A t = AC J™ F(| + ()sech 2 (y)dy « 2ACF((), 

C t = A 2 J™ F(| + CXsech 2 ^) - sech 4 (y))dy « ^ 2 F(C), (9) 

Q =2C + -J F(^ + ()ysech 2 (y)dy « 2C(1 + 0(— )). 
From these equations it follows that 00) = 1) 



a = 4,7(0, c fm = y&pf(CML). (10) 

In the next section we shall compare these predictions with direct numerical integrations 
of Eq.©. 

3. Numerical results 

To check these results we have numerically integrated Eq. (J2J) with the function 7(2) 
given by 



7(t) = 2 + A s 



Z 7T Z 



(11) 



This function models, for small s, an adiabatic change of the scattering length, while 
for large s it reduces to the step function of amplitude A s centered at I//2 (to have 
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Figure 2. Panel a) Bright and dark soliton width versus step amplitude A s , for 
parameter values s = 0.1, 7/ = 100, and soliton velocity v s — 0. The triangles 
and squares denote numerical values in absence of trap potential for a = 1 (bright 
case) and a = — 1 (dark case), respectively. The open circles denote numerical 
values for the bright case in presence of a trap potential of frequency uj 2 = 0.01 
(the soliton is taken at rest in the bottom of the trap). The continuous and dashed 
curves are the analytical predictions in absence of trap potential for, respectively, 
bright and dark soliton case. Panel b) The squared amplitude of a bright soliton 
vs center of mass position for different values of the amplitude A s of a kink-like 
spatial inhomogeneity centered at x=30, given by Eq. lllfl . The curves, from 
bottom to top, refer to A s = 2.5, 5.0, 7.5, 10.0, respectively. The other parameters 
are s = 0.2, Tf = 60, a = 1,Aq = 1. The soliton is initially at rest, placed at position 
Xini — 12.5. In the inset we show the soliton final velocity as a function of A s . The open 
dots are numerical values while the continuous curve is obtained from Eqs. I|10I10|) as 

Vfin = f^^-jl- 

the function starting from 2, the quantity 5 = 7(0) — 2 must be subtracted). In the 
following we shall use this function to model also spatial variations of the scattering 
length. In Fig. ^we depict the width of a bright (panel a) and a dark soliton (panel b) 
as a function of time for different values of A s , in the case u = 0. From this figure we see 
that the oscillations of the soliton width after the transition, are more pronounced for 
dark solitons than for bright ones. This is a consequence of the fact that in the dark case 
the background is also excited by the perturbation and contribute for a relevant part 
to the dynamics. In Fig. we show the average final value of the width of bright and 
dark solitons as a function of A s . It is evident that, for the same absolute variation of 
the scattering length, bright solitons are more compressed than dark ones, the difference 
being just a factor proportional to l/y^y, as expected from our analysis. The agreement 
with the theoretical predictions is indeed rather good especially for higher values of 
A s . Notice that in this figure the influence of a parabolic trap on the compression 
phenomenon is also shown for the case of a bright soliton initially at rest in the bottom 
of the trap (open circles). We see that the results are almost identical to those in absence 
of trap (this is quite expected, since the pulse is localized in the bottom of a low intensity 
trap). The panel (b) of the same figure shows the amplitude of a bright soliton as a 
function of the position of the center of mass for the case of a space dependent variation 
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Figure 3. Time evolution of bright and dark solitons in absence and in presence 
of parabolic trap. Panels (a,b): Bright soliton dynamics (a = 1) in presence of 
parabolic potential of frequency uj 2 = 0.04 subject to (a) an adiabatic change 7(f) 
of the scattering length with A s — 25, s ~ 0.5, and Tf — 50; (b) a step like variation 
7(t) with s = 10 4 . Other parameters are fixed as in panel (a). Panels (c,d): Dark 
soliton dynamics a = — 1 in absence (c) and in presence (d) of a parabolic trap of 
intensity uj 2 = 0.8. In both cases the initial dark profiles are at rest and subject to 
a 7(i) variation with parameters A s = 10, s = 0.5, and Tf — 200. Notice that the 
parabolic potential is also shown at t = to display the scale in the z-direction. 



of the scattering length in absence of the parabolic trap. The soliton, is initially at rest, 
is sucked into the higher scattering length region, and reaches a constant velocity after 
passing the inhomogeneity. Also in this case we find an excellent agreement with our 
analysis (see the inset of the figure). 

The effect of the parabolic trap on pulse compression is further investigated in 
Fig. El In panels (a) and (b) of this figure, we show the evolution of a bright soliton 
oscillating in the trap and subject, respectively, to adiabatic and abrupt variation of the 
scattering length. We see that in contrast with the adiabatic case, which always leads 
to compression, abrupt variations of the scattering length may induce the splitting of 
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the soliton into sub-solitons (notice in panel (b), that three solitons, indicated by the 
arrows, are formed at time t = 54). The splitting phenomenon is similar to the one 
observed in absence of parabolic potential 4J. The number of sub-solitons created can 
also be determined in analogy with this last case (we omit details for brevity). Panels (c) 
and (d) of Fig. H2 show the time evolution of a dark soliton in absence and in presence 
of a parabolic trap. We remark that, while in the bright case the compression of the 
soliton is obviously accompanied by an increase of its amplitude, in the dark case the 
amplitude cannot decrease below zero, so that the compression must be combined with 
a decrease of the background level (the matter from the background move inside the 
hole so to reduce its size). This is particularly evident in panel c for the case uj = 0. 
We remark that this case, although physically non realisable for BEC, it is of interest 
for nonlinear optics, and shows the role played by the confinement potential in the case 
of BEC with positive scattering lengths (it should be compared with panel(d) discussed 
below) . 

In analogy with bright solitons, one could expect that an abrupt change in time 
of the scattering length may split a dark soliton ground state into dark and grey sub- 
components if the amplitude of the perturbation is big enough. This is indeed what 
we numerically observed. For the case uj = 0, using the Inverse Scattering Transform, 
one can predict that out of an initial profile of the form ip(x,t = 0) = &tanh(or), one 
black soliton and a number iVo of pairs of grey solitons (with Nq — n — 1, n being the 
integer part of the ratio b/c) [T^HH)], should be generated. By denoting with 7x and 
72 the initial and final values of 7, one can write the initial condition for the equation: 
ivt + v xx — \v\ 2 v = 0, as v(0) = Woa/72 tanh(-u \/27 1 x), with v = ip^/y^, and the ratio 
b/c equal to A/72/271. We checked that this analysis indeed gives correct predictions of 
splitting in absence of the parabolic trap. The influence of the parabolic trap on a dark 
soliton initially at the rest in its bottom, is reported in panel (d) of Fig. El Although 
the phenomenon resemble the one observed in absence of the trap (notice the decrease 
of the background after the transition), the solution is more complicated due to the 
boundary effects introduced by the trap. In particular we see that the compression of 
the hole in the center is less effective than for uj = 0, and the decrease of the background 
is accompanied by an expansion of the condensate. This is a consequence of the non 
constancy of the background (it obviously must go to zero at large distances) which 
permits the matter to escape up in the trap (the pressure on the hole is thus partially 
lost in the expansion). The same phenomenon exists also for dark solitons initially 
oscillating in the trap. In analogy with the case uj = 0, we also found that abrupt 
variations in presence of parabolic potential induce dark soliton splitting. 

In Fig. 4 we show the averaged amplitude of the soliton (i.e. the time average of 
the modulo squared at the center of the profile) versus the parameter s controlling the 
rapidity of change of the function j(t). From this figure we see that a crossover from 
adiabatic to fast variations behavior occurs at s ~ 1.2. Notice that this crossover value 
of s corresponds to a variation function j(t) which is already quite rapid (see top inset of 
the figure). The adiabatic region is characterized by the fact that the final profiles of the 
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Figure 4. Time averaged amplitude of the soliton versus the parameter s of the 
scattering length function 7(f). The amplitude and the time of the variation are fixed 
as A s = 10 and Tf = 100. Notice the crossover from adiabatic to abrupt behavior 
occurring at s as 1.2 (dotted lines are drown to guide the eyes). The inset at the 
top shows the shape of the function j(t) at the crossover point (point b). The inset 
at the bottom shows the initial and final profiles of the soliton at the points a, 6, c, 
corresponding, respectively, to the values, s — 0.5, s = 1.2, s = 3.6. Soliton profiles 
corresponding to points a and c have been shifted by ±7.5 along the x axis to avoid 
overlapping with the ones corresponding to point b. 



compressed solitons are almost independent on s (they depend only on the amplitude 
of the scattering length variation). This can be seen from the bottom inset of Fig. 0]in 
which the initial (t = 0) and final (t = 100) soliton profiles corresponding to points a 
and b in the adiabatic region, are reported. Very little oscillations of the soliton profile 
are found in these cases. On the contrary, the non-adiabatic (or fast oscillation) region 
is characterized by the presence of large oscillations of the soliton profile, induced by 
the rapid variation of the scattering length. This is illustrated in point c of the lower 
inset of Fig. EJ in which we depicted, besides the initial profile at t — 0, two snapshots 
of the soliton (max and min of the profile) to evidenciate its oscillations at the final 
time Tf = 100. These oscillations increase by increasing s and eventually lead to the 
splitting of the soliton as illustrated in panel b of Fig. El (the splitting, however, depends 
also on the amplitude of the scattering length variation, and usually occurs much above 
the interface between adiabatic and fast oscillation behavior). From this analysis we 
conclude that the adiabatic compression of matter wave solitons is possible in a wide 
range of the parameter s (i.e. < s ~ 1) and is compatible with relatively fast variations 
of the scattering length, as one can see from the top inset of Fig. 4. 
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4. Conclusions 

In conclusion, we have shown that adiabatic changes of the scattering length, both in 
space and time, can be effectively used to control parameters of BEC matter waves 
of soliton type. In particular we showed the effect of pulse compression both on dark 
and bright solitons in presence and absence of a confining parabolic potential. The 
influence of the smoothness of the transition was also considered. We found that 
deviations from adiabatic predictions occur only for very large values of s, i.e. for 
almost step-like variations. This implies that the results of this paper are valid for 
generic smooth variations of the scattering length and therefore are expected to be 
found in real BEC matter waves experiments. The possibility of sharpening bright and 
dark solitons by means of smooth variations of the scattering length, beside providing 
another confirmation of the soliton nature of the BEC ground state, may be very useful 
for concrete applications. We hope that experiments in this direction will be soon 
performed. 
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